A finite element analysis of a silicon based double quantum dot structure 
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We present finite-element solutions of the Laplace equation for the silicon-based trench-isolated 
double quantum-dot and the capacitively-coupled single-electron transistor device architecture. This 
system is a candidate for charge and spin-based quantum computation in the sòlid state, as demon- 
strated by recent coherent-charge oscillation experiments. Our key findings demonstrate control 
of the elèctric potential and elèctric field in the vicinity of the double quantum-dot by the elèc- 
tric potential applied to the in-plane gates. This constitutes a useful theoretical analysis of the 
silicon-based architecture for quantum information processing applications. 

PACS numbers: 85.35.Be, 03.67.Lx 



Recent experiments conducted on trench-isolated dou- 
ble quantum-dot (IDQD) structures have successfully 
demonstrated detection of single-electron polarization, 1 
and coherent-charge oscillation>2 This highlights the pos- 
sibility of constructing charge-based quantum computer 
circuits in Si, with coherence times of the order 100 nsià4 
The architecture for a single qubit device is a com- 
plex, three-dimensional structure consisting of a single- 
electron transistor (SET), an IDQD, and gate elèctrodes. 
This makes it difficult to determine theoretically the sys- 
tem evolution by means of a complete and self-consistent 
Schròdinger-Poisson analysis. This is particularly the 
case if the analysis were to fully take into account the 
device geometry and all interactions while performing 
quantum manipulations within the coherence time of the 
qubit è 

In this paper, we present a significant contribution 
to such an analysis by a finite-element solution of the 
Laplace equation with the three-dimensional device ge- 
ometry and material composition taken into account. 
The aim of this work is to demonstrate the electrostatic 
effect on the IDQD structure when voltages are applied 
to the in-plane control gates of the device. 

Figuresnia) and^b) show device schematics in the x-y 
and cc-z-planes, respectively. The trench isolation, illus- 
trated in Fig.^b), is formed by high-resolution electron- 
beam lithography and reactive-ion etching. Each trench 
is approximately 150 nm deep and runs into the buried- 
oxide (BOX) layer of the silicon-on-insulator (SOI) wafer. 
The active regions of the device elements are P doped Si, 
which are electrically isolated from other device elements, 
as seen in Fig.^ In this work, we use rectangular approx- 
imations to the device elements and the etched profilcs 
to simplify the analysis & 

The small dimensions of the quantum dots and the 
20 nm constriction between them, which is fully depleted 
and acts as a tunable inter-dot tunnel barrier, result 
in a significant double-well typc confinement potential 
for the electrons that occupy discrete quantum states 
on each quantum dot£& The tunable inter-dot coupling 
causes the wave functions in the two dots to overlap and 
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FIG. 1: A schematic representation of the Si IDQD device 
used in the numerical simulations. Rectangular approxima- 
tions are made to the device elements, (a) A view of the x-y 
plane at the level of the dashed line in (b). (b) A view of the 
x-z plane at the level of the dashed line in (a). The trench- like 
structures are connected to the in-plane metal gates through 
a column of air. Each trench extends about 150 nm above the 
oxide base. 



hybridize so that they may be thought of as pseudo- 
molecular states of an artificial two atom molecule. The 
IDQD is an electrically-isolated component that is cou- 
pled only capacitively to the rest of the circuit, including 
the SET for read-out, and the in-plane control gates (Gl 
to G3) for manipulation. Voltages applied to the gates 
Gl to G3 are used to tune the elèctric field in the vicin- 
ity of the IDQD, and thus, the confinement potential 
asymmetry and inter-dot tunnelling. Hence, an electron 
initially localized on one quantum dot may be allowed to 
tunnel to the opposite quantum dot by such manipula- 
tion. 

An elèctric field is induced on the SET as a result 
of this polarization process. This modulates the chem- 
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ical potential of the SET, and, therefore, the conduc- 
tance through the source and drain leads under a finite 
SET bias condition. To maximize the change in conduc- 
tance, the SET is initially tuncd to the charge-sensitivc 
regime by an appropriate voltage bias at G4, and under a 
small source-drain bias to ensure operation in the linear 
transport regime. Such manipulation of the device has 
only been shown through experiment so far. Therefore, 
a thorough theoretical analysis is necessary to comple- 
ment the recent experimental findings, and build a more 
comprehensive understanding of the physical mechanisms 
involved. 

Analytic methods exist for calculating the electro- 
static potential in two-dimcnsional clcctron gases gen- 
erated by patterned surface gates on GaAs/AlGaAs 
heterostructuresAiS While these analytic methods yield 
useful results for such devices, they are unsuitable for 
trench-isolated Si structures, where the geometry is much 
more sophisticated. Therefore, numerical methods have 
to be implemented. The finite-element method is a well- 
suited means for simulation of geometrically-complicated 
domains^iíià and is commonly used to solve Poisson-type 
equationsii2ài4 

In order to determine the elèctric field throughout the 
modelled device regions, the numerical solution to the 
Laplace equation in thrcc-dimcnsions is performed: 



V • [e(x)V0(x)] 
<Xx) 

dn 



= 0, x G fi e M 3 , 
= Di, xe dü Di , 

= q, xe dü N , 



(1) 

(2) 
(3) 



where </>( x ) is the electrostatic potential and e(x) is the 
material dielectric parameter. The dielectric parameter 
varies discontinuously on moving through the different 
materials - air (e = 8.85 • 1(T 12 Fm" 1 ), Si (11.0e ) and 
SÍO2 (4.5eo). We apply Dirichlet boundary conditions 
to the surfaces of the metal gates and to the groundcd 
base of the device, and the Neumann boundary condition 
to the exposed surfaces. (We apply Neumann boundary 
conditions with q = 0, but for generality, we include in 
our discussion the possibility of non-zero q). 

The finite clement solution of the Laplace equation is 
well covered in the literatureiüüüà The bàsic idea of 
the finite clement method is to approximate the unknown 
fields, for example 4> m the Laplace equation above, by 
<f> which is a combination of linearly independent basis 
functions Na 
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where M is the number of basis functions and 4>j are 
the expansion coeficients to be determined. In the finite 
element method, the computational domain Çl is divided 
into a number of elements, and the Nj are chosen to be 



piecewise polynomials such that they are non-zero only in 
a 'few' adjacent elementsiüiíà The method then requires 
the substitution of <f> into the Laplace equation, and the 
residual R = V • [eV0] to be orthogonal to the space 
spanned by a linearly independent set {Wi, ■ ■ ■ , W n }. In 
our calculations, we implcmcnt the orthogonality through 
the weighted residual statement: 



RWidü = 0, 
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and Galerkin's method i.e. set Wj equal to the basis 
functions Ni to obtain 
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Using integration by parts, we reduce the order of 
derivatives in Eq. JHJ: 



d<j)dN t d^dN d( ^ dN i díï 
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Ne^dT = . (7) 

The weighted residual method and integration by parts 
leads to a natural mechanism for the incorporation of 
derivative boundary condition given by Eq. @ for the 
Laplace operator V • [eV</>]. Hence, we obtain after ex- 
panding the approximation for <f> 
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where summation is implied over repeated indices. The 
problem has now been reduced to one of matrix inversion; 



K<S> = F, 

where the 'stiffness' matrix K is given by 



(9) 



3N t dN 3 dNj dNj dNdN 
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içi dx dx dy dy 
and the RHS vector is given by 
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and $ is simply the vector of unknowns <f>j. Dirich- 
let boundary conditions are implemented by forcing pre- 
scribed vàlues of éj. 
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In the formulation and solution of Eq. (jÜJ) we have 
chosen linear basis functions corresponding to 8-noded 
brick elements. This method therefore has a convergence 
rate for error of 2.0H 

We have employed Gaussian quadrature in three- 
dimensions for the volume integrals and two-dimensions 
for surface integrals, and a conjugate gradient method 
with Incomplete Lower and Upper (ILU) factorization 
preconditioning to solve the resulting linear system of 
equations. Sufhcient resolution was obtained by a mesh 
with 107 x 77 x 9 nodes in the x, y and z directions, 
respectively. We implemented adaptive mesh refinement 
along z axis to improve accuracy in the vicinity of the 
active region. 
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FIG. 2: Two-dimensional slices taken from the full three- 
dimensional solution of the Laplace equation. (a) A slice from 
the x-y plane at z — 420 nm, through the center of the gates 
and IDQD. (b) A slice from the x-z plane at y — 190 nm, 
through the center of gate G2. The gate G4, is set to —4.8 V. 
The gates Gl, G2 and G3 are set to 1 V, 2 V and -1 V 
respectively. The position of the IDQD is outlined. 

Figures Efa) and|2h) show cross-sectional slices along 
orthogonal planes of the full three-dimensional solution 
for the simulated elèctric potential. The simulation was 
performed with the following parameters: the gate po- 
tentials of Gl, G2 and G3 are set to +1 V, 2 V and -1 V, 
respectively; G4 is set to —4.8 V. The voltage chosen for 
G4 is approximately equal to that used for this gate in 
the experimental demonstrations in Ref. ij where single 
electron polarization of the IDQD was obtained for such 
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FIG 3: Cross-sectional curves through the IDQD with x = 
320 nm and z = 420 nm (the dots of each curve correspond 
to nodal points of the computational grid). The gate G4 is 
set to -4.8 V. (a) The electrostatic potential of the voltage 
on gate G2 is changed from -5 V to +5 V. (b) The voltage 
applied to gate Gl is varied from V to -5 V. The voltage 
applied to gate G3 is set to the negative of that applied to 
gate Gl (in order to maximize the elèctric field). 



a device. The choices of Gl, G2 and G3, are also similar 
to those in the experiments but for this simulation, the 
exact vàlues are chosen so that three-dimensional illus- 
trations are as clear as possible. 

Figures |5fa) and|2fb) clear ly demonstrate the effect of 
the applied gate voltages on the potential landscape of 
the IDQD and the device as a whole; the result of apply- 
ing a voltage on the in-plane gates is that a significant 
fraction of the applied voltage is induced on the IDQD, 
despite the etched trench gap. The abrupt change in the 
cffcctivc pcrmittivity from the metallic gates to the voids, 
from the voids to the SÍO2, and from Si to SÍO2, causes 
somc dcfinition of the gates and the IDQD in the plots. 
The difference between the relative permittivity of air, 
Si and SÍO2 leads to a potential gradient, such that the 
absolute value of the potential is prone to vanish more 
rapidly in air, compared with Si and SÍO2. However, 
Fig. Efb) clearly demonstrates that for this particular 
pillar height, which matches the device used in exper- 
iment, the potential at the IDQD is due mainly to the 
elèctric field vectors that are on a direct path through the 
trench isolation, and not the underlying substrate. This 
is consistent with experimental observations and is the 
preferred mechanism of device operation, since it is rela- 
tively easier in design and theoretical analysis, compared 
with the case where the majority of the elèctric field is 
through the semiconductor base and the field lines arrive 
at the IDQD from several different paths. 

For a more quantitative analysis, we determine the 
elèctric potential along the active region of the IDQD 
as a function of the applied gate voltages. This is shown 
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in Figs.|2ía) and[3{b), where difïerent gates are used to 
apply the in-plane elèctric field. Figure^a) shows that 
the effect of varying the voltage applied to gate G2, from 
-5 V to +5 V, is to induce a voltage at the IDQD from 
-2 V to +0.8 V, respectively. Note that G4 had lowered 
the overall potential by approximately 0.55 V in this case. 

The data in Fig. EJa) shows a maximum change of 
0.3 V of the electrostatic potential at the IDQD, when 
the G2 gate voltage is raised or lowered by 1 V. This 
field coupling factor of ^30 % is approximately one or- 
der of magnitude greater than what was observed in 
experimenti However, the measured quantity in experi- 
ments is the SET current, and the IDQD coupling terms 
are inferred from such measurements. The task of cal- 
culating such coemcicnts exactly as measured in experi- 
ment is beyond the scope of a purely electrostatic model, 
since, with the SET present, the global system that must 
be treated consists of intcracting sub-systems of quantum 
mechanically bound electrons. Therefore, we project that 
a self-consistent Schròdinger-Poisson analysis of the sys- 
tem would yield results for the coupling coeficients that 
are closer to actual vàlues observed in experiment. 

Our results are also consistent with the experimental 
demonstrations of Ref. Q where the voltage on gate G2 
is swept continuously from —5 V to +5 V, with G4 set 
to ^—4.8 V, in order to demonstrate conductance res- 
onances of the SET currents but also resonances due to 
the single electron polarization of the IDQD. (The device 
in Ref. [lj had a slight asymmetry in the alignment of the 
IDQD relative to G2, hence the ability of G2 to polarize 
the IDQD.) 

The abrupt changes in the potential at 115 nm and 



265 nm are due to the change of relative permittivity at 
the air-Si02 interface. The difference between the poten- 
tial gradient in the air and in the semiconductor regions 
is morc evident in thcsc figures. Figure IHJb) shows the 
effect of applying voltages of opposite sign to gates Gl 
and G3. This results in a potential gradient across the 
IDQD, which has a maximum value of ~0.007 Vnm -1 in 
these simulations. This clearly demonstrates an effective 
mechanism for externally tuning the infernal potential 
asymmetry of the IDQD electrònic states. 

In the experimental demonstrations of Ref. |2], a volt- 
age bias is pulsed across in-plane metallic gates, which 
are placed perpendicularly to an IDQD as in our case, 
and this was shown to result in the coherent oscillation 
of a single electron charge present in an IDQD. This is 
again consistent with our results which suggests a strong 
elèctric field is induced at the IDQD due to the elèctric 
field at the in-plane gates. 

In summary, we have successfully demonstrated, by 
means of finite-element solutions to the Laplace equation, 
that the elèctric potential and potential gradient across 
the confining region of the IDQD in trench-isolated Si 
devices may be manipulated effectively by the voltages 
applied to capacitively-coupled in-plane gates. Our cal- 
culations show good correlation with recent experimental 
demonstrations, where the IDQD electron states are ma- 
nipulated by such methods. 
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